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Abstract 

We examine the shapes and energies of 5- and 7-fold dischnations in low-temperature 
hexatic membranes. These defects buckle at different values of the ratio of the bending 
rigidity, k, to the hexatic stiffness constant, Ka, suggesting two distinct Kosterlitz-Thouless 
defect proliferation temperatures. Seven-fold disclinations are studied in detail numerically 
for arbitrary k/Ka- We argue that thermal fluctuations always drive k/Ka into an "un- 
buckled" regime at long wavelengths, so that disclinations should, in fact, proliferate at the 
same critical temperature. We show analytically that both types of defects have power law 
shapes with continuously variable exponents in the "unbuckled" regime. Thermal fluctu- 
ations then lock in specific power laws at long wavelengths, which we calculate for 5- and 
7-fold defects at low temperatures. 
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1 Introduction 



Topological defects, such as dislocations and disclinations, behave differently in membranes 
and monolayers. Monolayers are films, typically with a triangular lattice in their low- 
temperature crystalline phase, strongly confined to a plane by, for example, surface ten- 
sion. In monolayers, point-like topological imperfections are the ingredients of a detailed 
two-stage melting theory, which predicts that the usual latent heat associated with a 
first order melting transition can be spread out over an intermediate hexatic phase char- 
acterized by long range bond orientational order and short range positional order. 0, |^ In 
contrast to monolayers, membranes can easily buckle out of the two-dimensional plane. 
This buckling lowers defect energies. The standard example of a membrane is an extended 
lipid bilayer surface in water solution. 0] Dislocation energies in membranes are finite, in 
contrast to a logarithmic divergence with system size for monolayers, which leads to melt- 
ing of crystalline membranes at any nonzero temperature for entropic reasons, Q The 
generic low temperature phase for membranes at large enough length scales is inevitably a 
hexatic liquid with long range bond orientational order and a vanishing shear modulus. 

Disclinations can also lower their energy by buckling. The energy of a single discli- 
nation in an otherwise crystalline monolayer diverges with the square of the system size. 
Buckling in a membrane leads to energies of plus and minus disclinations that diverge only 
logarithmically with size.[^] Interestingly, the coefficients of these logarithmic divergences 
are different: the energy of 5-fold (plus) disclinations is about a factor of 2 lower than that 
of a 7- fold (minus) ones.0 This asymmetry differs from the behavior of defects in most 
monolayer materials that exhibit the Kosterlitz-Thouless defect unbinding transition. 
Plus and minus vortices in superfiuid helium films, for example, must have exactly the 
same core energies and logarithmic divergences with system size by time reversal invari- 
ance. Equality of energies also holds for dislocations with equal and opposite Burgers 
vectors in crystalline monolayers. Plus and minus disclinations in hexatic monolayers have 
different core energies, due to the different local environments of the 5- and 7-fold defects, 
but the coefficients of their logarithmically diverging energies are identical. []I|] Identical 
coefficients for the logarithmically diverging ± defect pairs ensures that the elementary 
Kosterlitz-Thouless energy-entropy balance leads to the same proliferation temperature for 
each type of defect. This famous argument predicts that 5- or 7-fold defects prolifer- 
ate whenever the free energy to create a disclination, F^lR) = E^,{R) — 2k-BTln{R/ao) or 
F-j{R) = Ej{R)—2k-QTln{R/aQ), becomes negative. Here E^{R) and Ei{R) are disclination 
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energies as a function of the system size R, and ao is a microscopic cutoff. 

The disparate dischnation energies in buckled crystalhne membranes suggests that these 
energies may also differ in hexatic membranes, as recently emphasized by one of us.[^ 
The energy of 5-fold disclinations in hexatic membranes has been studied by Seung,P| 



Guitter, and Guitter and Kardar. |Tl| The results depend on two dimensionless parame- 
ters, K/k-gT and KA/k^T, where k and Kj^ are the membrane bending rigidity and hexatic 
stiffness constant, respectively. When k/K^ ^ 1, the membrane remains asymptotically 
flat in the presence of both 5- and 7-fold defects, and the energies diverge logarithmically 
with equal coefficients. When k/Ka < 11/72, however, the 5- fold disclination buckles and 
reduces its in-plane bond angle energy at the expense of bending energy. The resulting 
energy still diverges logarithmically, but the coefficient is reduced by buckling. The locus 
of Kosterlitz-Thouless transition temperatures for 5-fold defects in membranes when buck- 
ling is included is shown by the solid curve in Fig. 1. These disclinations assume a conical 
shape when k/Ka < 11/72 and, as shown in Sec. 3, exhibit a weaker power law deviation 
from flatness when k/Ka > 11/72, with continuously variable exponents. 

In this paper, we study 7-fold defects and determine the curve for their proliferation, 
shown as the dashed curve in Fig. 1. Because buckled 7- fold defects do not have the 
azimuthal symmetry of 5-fold defects, some numerical work is necessary. As anticipated 
in [|], the 7-fold hexatic disclination energy exceeds its 5-fold counterpart for k/Ka -C 1, 
leading to a larger Kosterlitz-Thouless critical temperature. We find that 7-fold disclina- 
tions buckle to reduce the coefficient of their logarithmically diverging energy whenever 
k/Ka < 13/216, while they remain asymptotically fiat when h/Ka exceeds this value. 
There is again a power law deviation from flatness for h/Ka > 13/216. As is evident from 
Fig. 1, naive application of the Kosterlitz-Thouless criterion to hexatic membranes leads 
to two distinct defect proliferation temperatures over a significant range of parameters. 

Our results provide useful information about deformations of hexatic membranes near 
defects at relatively low temperatures. However, understanding the behavior at very large 
distances from the defect cores requires that we take into account the nonlinear renor- 
malization of n and Ka by thermal fluctuations. These effects were first studied in a 
perturbative, low-temperature expansion by F. David et al.,^^ who found that the hex- 
atic phase of membranes is controlled by a line of fixed points in the {k-QT / KAik-QT / k) 
plane with slope k* /K\ = 1/4. The parameters k and Ka are driven toward the fixed line 
by thermal fluctuations even if they initially lie in one of the buckled regimes. This line 
and the associated nearby renormalization group flows are indicated by the heavy curve in 
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Fig. 1. The unstable Kosterlitz-Thouless fixed line for hexatic monolayers (k oo) is also 
shown. Because the stable fixed line has a slope far below the critical slopes for buckling 
of 5- or 7-fold disclinations, we conclude that these defects probably have symmetrically 
diverging energies when thermal fluctuations are taken into account. Provided that the line 



of stable fixed points, which is only known perturbatively at low temperatures, |T^ does 
not bend significantly before piercing the vertical part of the Kosterlitz-Thouless instabil- 
ity curve, both defects will remain unbuckled and proliferate at the same point. Although 
the coefficients of the logarithmic divergences in 5- and 7-fold disclination energies are the 
same, these defects have interesting power law shapes in the "unbuckled" regime. 

Park and Lubensky have recently incorporated fluctuations in the internal metric of the 
membrane into the work of David et a/.|]T3 Inclusion of these fluctuations appears to be 
necessary to account for local in-plane shear and compression modes, whose primary effect 
is to dilate and reshuffle the nearest-neighbor bond connectivity of the atomic or molecular 
constituents. Averaging over these modes leads to an effective coupling K'^ which should 
replace the hexatic stiffness in the model used here. The bending rigidity is unchanged. 
The remaining renormalization of k and K'^ by thermal undulations, however, is identical 
to that found by David et al. Thus, our overall conclusions are unchanged, provided we 
use K and the renormalized coupling K'^^ in our results and in Fig. 1. In agreement with the 
results of the study of disclination statistical mechanics by Park and Lubensky, ||13|| we still 
expect a symmetrical ± disclination unbinding transition with unbuckled defects in hexatic 
membranes. Disclinations cannot unbind separately, in contrast to the predictions of the 
"naive" Kosterlitz-Thouless argument, provided the thermal renormalization of the ratio 
i^/K'^A to 1/4 at long wavelengths remains intact out to the transition temperature. |I5| 

In independent work. Park and Lubensky have also studied the buckling of 7-fold 
disclinations. [jl4| Their result for the critical ratio n/ Ka and the behavior near the buckling 
transition are in agreement with the results presented here. Our results are more reliable 
in the limit n/K^ ^ 1, however. 

It is worth noting that similar issues involving disclination asymmetry arise even for 
monolayers, when two-dimensional nematic order is present. The topologically stable de- 
fects are now ±1/2 disclinations, and the ordered state is described by both bend and splay 



elastic constants Ki and /^s.lTBI When Ki ^ K^, the energies of ± disclinations again di- 
verge logarithmically with different coefficients. Here, a renormalization group analysis of 
nonlinear effects due to thermal fluctuations shows that Ki and are driven to equality 
at long wavelengths, so that the "one-Frank-constant approximation" becomes asymptoti- 
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cally correct at long wavelengths, [p!^] The disclinations energies are equal in this limit, and 
one again expects a symmetrical defect unbinding transition. [|l^ 

We should stress that even though disclination energies are asymptotically symmetrical 
in hexatic membranes, buckling will still occur locally for appropriate parameter values. 
Buckling will persist out to length scales such that the renormalized value of the ratio 
k/K'^^ is in the unbuckled regime. An intrinsic microscopic asymmetry, moreover, can still 
arise in the liquid, after the disclinations unbind. As emphasized in Ref. 0, a net excess of 
disclinations should have important consequences in liquid membranes with free boundary 
conditions or a topology that can change on experimental time scales. Exactly how the ± 
disclination populations become identical as one approaches the liquid-to-hexatic transition 
from the liquid phase is an interesting topic for future research. 

In Sec. 2 we discuss how membrane buckling can screen disclinations. In Sec. 3 we review 
analytical calculations of the energy of a 5-fold disclination. We use approximate theory 
and exact numerical calculations to calculate the energy and shape of a 7-fold disclination 
in Sec. 4. We review the important effects of thermal fluctuations, including how these fix 
the asymptotic defect shapes in the unbuckled regime, in Sec. 5. 



2 Curved Hexatic Membranes with Defects 

The Hamiltonian for a flexible, hexatic membrane is given in the limit of vanishing surface 
tension by P, |19| 



Hj, = kJ d'SK . (1) 

All these integrals are over the surface of the membrane. For the case of a membrane with 
free boundary conditions there should also be line tension and geodesic curvature terms. 
We neglect these terms. Here is the inverse metric tensor, H is the mean curvature, K 
is the Gaussian curvature, and the gauge field is defined by i^' = Y^DiQj. Upon defining 
g = det{gij), we have = eijg^^^"^ with en = £22 = and eu = —£21 = 1- The surface 
area element is given by (PS = cPag^^"^. 

We are interested in very flexible membranes, as opposed to monolayers, and so we will 
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neglect the possible surface tension term of the form Hj. = r J d?S. We will also ignore the 
Gaussian curvature term if^, which is a perfect derivative by the Gauss-Bonnet theorem. 
i^K is the standard bending energy term, and Hg is the contribution to the energy from 
fluctuations in the local bond order parameter. The bond order parameter is frustrated by 
the rotation of tangent vectors that occurs under parallel transport on a curved surfaces. 
The amount of frustration is given by the gauge field, 

To gain some physical understanding of a flexible hexatic membrane, we examine the 
ground states. In particular, we search for the low temperature geometries of 5- and 7-fold 
disclinations. After performing the minimization over 9, we find 

5H 

se{a) 

where 6q is the bond angle field that minimizes the energy. Upon defining 



= g-'/^[g''{^A-n^)] = 0, (2) 

0=00 



diOo -n^ = -fi^djX , (3) 

we find the derivative (2) is zero except when derivatives of x do not commute. We can 
show that disclinations cause the derivatives to fail to commute by applying the operator 
7'='5fcto Eq.gg 

D.D'x = K-^^'dudA 
= K{a)-s{a) 

= c(a) . (4) 

The disclination density is given by 

s{a) = Y.s,S{a-aM^^r'^' ■ (5) 

i 

Here the disclination strength is given by Sj = 7r/3 for a 5-fold disclination and Sj = — 7r/3 
for a 7- fold disclination. Given the form of Eq. (^, we can express x terms of a Green's 
function as 

X{a) = J d'S'G{a\a')c{a') , (6) 

where 

D,D'Gia\a') = V^G(a|a') = 5(a - a')gia')-'^' . (7) 
The Hamiltonian when the bond angle field is minimized is given by 

He = ^l (fSd'S'(fS"c{a')c{a") [diG{a\a')]g'^ [d,G{a\a")] . (8) 
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We see Eq. that the relevant quantity is not the dischnation density or the Gaussian 
curvature separately, but rather the difference, c(cr), between them. Consequently, the 
hexatic energy arising from a disclination can be reduced by a non-zero Gaussian curva- 
ture. This screening, of course, will cost the membrane in terms of bending energy. The 
competition between screening of the hexatic energy and bending energy determines the 
equilibrium shape of the membrane. 

The energy of a single, isolated disclination with "charge" s in a flat, circular membrane 
is given hj E = {s'^/Att) ln(i?/ao).[|ll Here R is the radius of the membrane and Qq is a mi- 
croscopic cutoff. For 5- and 7-fold disclinations, = Ei = (vrii'yi/SG) ln(i?/ao). Buckling 
of the membrane can reduce this energy. We describe the location of the membrane by 

X(r, 0) = (r sin 0, r cos 0, /(r, 0)) . (9) 

The diverging contribution to the energy comes from the large r region of the surface. The 
bending energy can diverge no more strongly than log(i?/ao) in a buckled ground state, 
since otherwise the energy would increase upon buckling. This bound implies that / grows 
at most linearly with r. If / grows less rapidly than r, then the Green's function defined by 
Eq. (1^ is given by drG ~ 1/ (27rr) as r ^ oo.p^ Furthermore, the Gauss-Bonnet theorem 
then implies that / (PSK = for membranes with a disk-like topology. I^T] From Eq. (|^), 
we see that the hexatic energy remains E ~ (tcKa/^Q) ln(i?/ao) and has not been reduced. 
For the logarithmic hexatic energy to be screened by buckling, therefore, the height must 
grow linearly with r: 

X(r, 0) = (rsin0, r COS0, r/i(0)) . (10) 
The Green's function that satisfies Eq. (0) is then given by 

G(r,0) = rMn[r(l + /i(0)2)i/2] , (11) 

with 

To evaluate the hexatic Hamiltonian, we need both the Gaussian and mean curvature. 
For the surface defined by Eq. (^, the Gaussian curvature is proportional to a delta func- 
tion: 

K{a) = a6{a)g-^'^ . (13) 
The coefficient, a, can be determined from 

a = ( d^SK 
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da'Vti , (14) 

where c is a contour bounding the surface. If ei and e2 are an orthonormal basis for vectors 
tangent to the surface, the gauge field is given by [p 



fij = ei ■ 9^62 . (15) 

From this equation we find, by taking ei and 62 to be basis vectors in polar coordinates 
and subracting the result for a flat surface, 



fii = 

(1 + /l2 + /i'2)l/2 



We, therefore, conclude that 

a = 2Tx-h. (17) 
We can now perform the integrals in Eq. (|^) to find 

= - l..(fl/«o) . (18) 

The mean curvature is given by 



(1 + |V/|2)V2 
{h + h"){l + h^) 



r(l + /l2 + /l'2)3/2 • 

The bending energy is, then, given by 



(19) 



The contribution to the bending energy associated with the singularity at r = will be 
absorbed into a core energy. The total energy of a hexatic membrane with a single, isolated 
disclination, excluding the core contribution, is given by 



ln(i?/ao) 26 2J ^ (1 + /i^ + /i'2)5/2 

The geometry of lowest energy is found by minimizing with respect to the function /i(0). 
Note that h depends on hi^cj)) though Eq. (0). 
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The route from the covariant Hamiltonian (|I]) to the tractable expression (^Tj) is com- 
phcated. For a nearly flat surface, a simplified Hamiltonian is often used H 

H=^jdh [d.e - + I / rfV(VV)^ , (22) 

where 

A = \e^Mdif){d^f)] , (23) 

and the derivatives are in flat space. The bond angle field, ^0(0") that minimizes this energy 
is given by 

9,(9A-^.) = 0. (24) 

To satisfy this equation, we define 

{dA - A,) = eijdjX ■ (25) 
Applying the operator eikdk to this equation, we find 

V\ = {dlf){dif) - {d^djr - s{r) , (26) 

with 

s(r) = ^s,5(r-r,) . (27) 

i 

For an isolated 5- or 7-fold disclination, s(r) = ±(7r/3)(5(r). The Hamiltonian now reduces 
to 

H = ^J rfVlVxr + f / rfV(VV)^ . (28) 

We can further find the height function, /, which minimizes this Hamiltonian. It satisfies 
a second nonlinear, hexatic "von Karmon equation" : 

^VV = idixMf) + {dlx){dlf) - 2{d^dyx){d.dJ) . (29) 

With the simple Hamiltonian (|2^), then, we have explicit partial differential equations 
that define the surface of minimal energy. For the covariant Hamiltonian, the differential 
equation arising from minimizing Eq. (^1]) is much more complex. 

3 The Energy of a 5-fold Disclination 

A 5-fold disclination can be screened by a surface with a positive Gaussian curvature. The 
natural surface to consider is a cone. 
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We first review the results of tlie approximate Hamiltonian (|22|). A cone defined by 
/(r) = ar solves Eqs. (|26|) and (|29|) , with x(r) = — (k/T^^) ln(r/ao).0 The coefficient is 
given by = 1/3 — 2k/Ka- For k/Ka < 1/6, the membrane buckles. The energy is given 
by 

f (vrK/3)(l - 3/t/KA)ln(i?/ao), /t/i^A < 1/6 
^5 ~ ■ (30) 

[ (7ri^A/36)ln(/?/ao), «:/ifA > 1/6 
We now review the results of the covariant Hamiltonian ([^).||T^ Equation ( pT]) fully 
specifies the energy, with h = a and b = 271(1 + a^)~^/^. We first note that Eq. ( p!7D can be 
derived from a geometrical argument. We consider capping off the cone with a small sphere 
of radius e, as in Fig. 2. The bending energy is unaffected by this small perturbation, since 
we are ignoring the contribution near r = 0. The Gaussian curvature is zero everywhere 
except on the sphere. On the sphere it is given by a = / <PSK = 2n J^^^^ du = 27r(l— cos ip). 



With tan?/) = a, we have a = 27r[l — (l+a^) ^/^], in agreement with Eq. ([T7|). Upon defining 



X 



'1 + a^) we have 



H 



(l/6-l + a;)2 ^ K 1 



X 



X 



X 



ttKa ln(i?/ao) 
Minimization of this equation leads to 

' tiKa {2 [(25/36 + k/Ka){1 - k/Ka)]^''' - 5/3} 
X ln(i?/ao), k/Ka < 11/12 



(31) 



(32) 



(7ri\:A/36)ln(i?/ao), k/Ka> 11/12 
We note that the limit Ka oo corresponds to the inextensional limit of a crystalline 
membrane. This energy has the correct limit (ll7rK/30) ln(i?/ao) as Ka —>■ oo, which 

corresponds to a crystalline membrane. |^ While the Hamiltonian (|2^) is often thought of 
as valid for small V/, we see that it does not exactly predict the buckling transition, where 
V/ is a small, nonzero constant. 

When k/Ka > 11/72, the above calculation shows that the height grows sublinearly 
with r. In fact, we now show that the height grows with a power that depends continuously 
on ti/ Ka- We assume that /(r, 0) = /(r). The hexatic and bending energies of Eq. (|l|) are 
then given by 



-kK, 



dr- 



f 



/2\l/2 



r 



1 



+ 



rf" 



(l + /'2)l/2 (l + //2)3/2 



(33) 
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We have set the short-range cutoff to oq = 1. As we shall see, /' — as r ^ oo, when 
i^/Ka > 11/72. Upon expanding Eq. (|33|) for small /', we find 

/oo 
dr 

Upon solving the equation 6F/6f{r) = 0, we find in the case of constant moduli 

/(r) ~ ar^^^ as r oo , (35) 

with 

y=[l-llKA/{72K)]'/^ . (36) 

We also minimize Eq. (|33| ) numerically. We express f'{r) on a grid at grid points = 
exp(lnrjnaxV''^) approximate the integral by a sum and derivatives by finite difference. 
We found convergence was achieved for n = 200, rmax = 100. Results are presented for 
n = 400,rinax = 200. The Polak-Ribiere conjugate gradient method was used to determine 
/'(rj).|2^ Figure 3 shows the height as a function of r for the specific case k/Ka = 1/4. The 
numerical results reproduce the asymptotic scaling of Eq. (P^j). The energy of this ground 
state is = {ttKa/^Q) InR — 0.0281'n'KA + E^, where E^ is a core energy contribution. If 
this core contribution is sufficiently large, Ec > 0.028l7ri^^, the surface of minimal energy 
would be flat, and the constant a in Eq. (|35|) would be zero. 



ir\2 



Ka ^ f^if' + rf") 
36r r 



72r 



(34) 



4 The Energy of a 7-fold Disclination 

A 7-fold disclination can be screened by a surface with a negative Gaussian curvature. There 
is no obvious natural surface to consider in this case. Using the approximate Hamiltonian 
(p2D, we can achieve an analytical answer, however. We let 

/(r, 0) = ar sin 20 
3k, 

X{r) = — ln(r/ao) . (37) 

Equation ( pUD is solved by this choice. Equation (|2^) is solved provided = 2/9 — 4^/ Ka- 
When k/Ka < 1/18, the membrane buckles. The energy is given by 

(7r/t)(l - 9k/Ka) HR/ao), k/Ka < 1/18 
E,^{ . (38) 

(7rir^/36)ln(i?/ao), k/Ka> 1/18 
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The covariant Hamiltonian ([1|) does not yield so easily to an analytical treatment. We 
can, however, expand Eq. (^) for small h to find the buckling transition: 

= = ^^{h + h") + K{h + 2h" + h"") + 0(/i3) . (39) 
on{(p) 72 

This equation predicts buckling for n/K^ < 13/216 with = a sin 20. Again the 

approximate Hamiltonian (^2]) predicts the transition value only approximately. 

We determine the surface that minimizes the energy for general values of i^/Ka by 
numerically identifying the function h{(j)) that minimizes Eq. (0). We express h{(f)) on a 
grid at grid points 0j = 27ri/n, again approximating the integral by a sum and derivatives 
by finite difference. We found convergence was achieved for n = 100. Results are presented 
for n = 200. The Polak-Ribiere conjugate gradient method was used to determine the 

Figure 4 presents the numerically-determined energies Hg, H,^, and H. The membrane 
buckles when h/Ka < 0.060, in good agreement with the exact value of 13/216 As the 
membrane becomes less stiff, the buckling is able to screen more and more of the hexatic 
energy: the hexatic energy goes from 

He ~ (vrKA/36) \n{R/ao) for k/Ka > 13/216, (40) 

to 

Hgr^O as k/Ka 0, (41) 
as expected. Similarly, the bending energy goes from 

H^r^O for k/Ka > 13/216, (42) 

to 

~ 2.27/tln(i?/ao) as k/Ka ^ . (43) 

The limit as k/Ka agrees with numerical calculations for inextensional crystalline 
membranes. 

Figure 5 shows the function h{(f)) for various values of the ratio k/Ka- As expected, 
the surface is flatter for stiffer surfaces. For very flexible membranes, the surface converges 
to a limiting shape. This limiting shape is very nearly proportional to sin 20, as shown in 
Figure 6. More generally, we can expand /(r, 0) in a Fourier series 

oo 

E /m(r)cos(m0) , (44) 

m=0 
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where /m(0) = 0. All odd terms must vanish for a two-fold symmetric saddle point con- 
figuration. In addition, / should change sign under a 7r/2 rotation, which eliminates the 
terms in Eq. (|^) with m = 0,4, 8, . . .. Such a symmetric saddle has an expansion of the 
form 

oo 

fir, 0) E U+2{r) cos [2(2p + 1)0] , (45) 

a conclusion also reached by Park and Lubensky.[P^]. We have checked numerically that 
the only non-zero Fourier components in Eq. (^) are indeed of the form m = 4p + 2, 
although the m = 2 term alone provides an excellent approximation. 

When k/Ka > 13/216, the height grows sublinearly with r. Just as for the 5-fold 
disclination, the height grows with a power that depends continuously on k/Ka- To see 
this, note first than when k/Ka > 13/216, V/ ^ as r oo. Upon expanding Eq. (|T|) 
for small V/, we find 

Ka 



^ / drd(j) 2/r + 13{drfy/r + ll(9^/)7^' 



12 idrid^fY /r 



+ -jdrd<Pr[dtf + {drf)lr + {dlf)/r'\ + Oif 
The solution of 6F/6f{r, 0) = is 



(46) 



fir) 



ar 



i-j/ 



sin 20 as r — > cx) , 



(47) 



with 



y 



720 + 13Ka/k - (331776 + 29952Ka/k + imK\l k 



2U/2' 



1/2 



/12. 



(4J 



As in the case of 5-fold disclinations, this sublinear decay leads to an additive, constant 
correction to the logarithmically diverging energy as i? ^ oo. 



5 Thermal Fluctuations 

We have so far ignored thermal fluctuations of the hexatic membrane. This assumption is 
valid only for the T — > limit. For finite temperatures, and for large membranes, thermal 
fluctuations will become important. 

A complete discussion of thermal effects is beyond the scope of this paper. We can, 
however, use the results of David et al. [O] and of Park and Lubensky |TB[ to estimate how 



the structure of disclinations in hexatic membranes is modified at finite temperatures. Park 
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and Lubensky argue that proper implementation of an ultraviolet cutoff to fluctuations in 
hexatic membranes leads to the replacement 

Ka ^ Kf = Ka- ^k^T{KA/t^f . (49) 

The bending rigidity k is unchanged. In the absence of a non-zero disclination density, the 
remaining renormalization equations for n and K'^ are those found originally by David et 
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dKf 
dl 



dK/ksT _ _ 3 / K 



cfT^ 
A 



dl 47r V Ak 



(50) 



The renormalization group flows induced by these equations are indicated schematically by 
the arrows in Fig. 1. 

We apply these results to the dilute limit of isolated disclinations discussed in Sees. 4 
and 5. The locus of disclination unbinding transitions, given by the criteria -F5(k, Ka) = 
and Fj{k,Ka) = discussed in the Introduction, are shown as the solid and dasked lines 
in Fig. 1. When thermal fluctuations are superimposed on the solutions of the T = 
extremal equations for disclinations in a membrane of size R, standard flnite size scaling 
arguments suggest that the couplings controlling the defect energies on this scale should 
be the running coupling constants K'a{1) and /?(/) obtained from Eq. ( pOD evaluated at 
I = ln(i?/ao). Effects of thermally-excited bound disclination pairs on an otherwise isolated 
defect could be included by adding a vortex fugacity to the set of recursion relations. We 
then expect that Eqs. (|33|) and should be replaced by expressions where Ka and n are 
replaced by the functions K'a{1 = Inr/oo) and k{1 = Inr/ao) appearing inside the integrals 
over r. Although K'^ does not renormalize at this order, a nontrivial renormalization could 
appear when higher order corrections in k^T/K and k-^T /K'a are included. 

If the basin of attraction of the locally stable flxed line in Fig. 1 includes the entire 
hexatic phase, K'A{r) and /t(r) will always be driven as r ^ cxd into the unbuckled regime 
for both 5- and 7-fold disclinations, since lim^^oo '^('")/-^a^('") = 1/4. The bending energy 
will then not contribute to the logarithmically diverging part of the energy, and both 
defects should unbind at the same point. Although "unbuckled" in this sense, the limited 
disclination shapes will be characterized by the power laws (35) and (0), with n/ Ka = 1/4. 
For 5-fold defects we flnd the asymptotic shape is given by ?/ = ?/+ = 0.6236, while for 7-fold 
defects we have y = y- = 0.8249. 
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Figure Captions 

Figure 1. The phase diagram for proliferation of isolated 5-fold (solid) and 7-fold 
(dashed) disclinations. Five- and seven-fold disclinations buckle above lines (not shown) 
extending from the origin to the tops of the vertical portions of the solid and dashed 
curves, respectively. For small k, or Ka (outside the curves), disclinations proliferate. 
Lines of renormalizcd effective rigidities are also indicated (bold). Renormalization group 
flows obtained by David et al. in a low temperature perturbation expansion away from 
the unstable fixed line at = oo to the stable line describing the crinked phase are 
indicated by the arrows. 

Figure 2. The cone /(r) = ar shown in projection capped by a small sphere used to 
calculate the integrated Gaussian curvature. The angle ip is given by tanip = a. 
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Figure 3. The height of an "unbuckled" membrane with a 5-fold disclination as a function 
of r for the case h/Ka = 1/4. There is a short distance cutoff so that /(r) is undefined 
for r < ao- 

Figure 4. The hexatic (short dashed), bending (long dashed) and total (solid) energies 
when r — s> oo as a function of Ka for a 7-fold disclination. An overall factor of ln(i?/ao) 
has been suppressed in each term. The defect is unbuckled for n/ Ka > 13/216 ~ 0.060. 

Figure 5. The surfaces /i(0) for a 7-fold defect above its unbuckling transition for the 
cases k/Ka = 0.06, 0.05, 0.03, 0.01, and 0.001. 

Figure 6. The surface /i(0) in the limit k/Ka — > (solid), which mimics the behavior in a 
crystaUine sohd, and the function 0.534 sin(20) (dashed) for a 7-fold defect. 
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